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We extend the notion of "quantum blob" studied in previous work 
to excited states of the generalized harmonic oscillator in n dimensions. 
This extension is made possible by Fermi's observation in 1930 that 
the state of a quantum system may be defined in two different (but 
equivalent) ways, namely by its wavefunction 'J or by a certain function 
^> ■ 5f on phase space canonically associated with VP. We study Fermi's 

5^ , function when "J is a Gaussian (generalized coherent state). A striking 

result is that we can use the Ekeland-Hofer symplectic capacities to 
characterize the Fermi functions of the excited states of the generalized 
harmonic oscillator, leading to new insight on the relationship between 
symplectic topology and quantum mechanics. 
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1 Introduction 



1.1 What We Want to Do 

We address the question whether it is possible to represent geometrically a 
function ip of the variables x = (x±,X2, •-, x n ). The problem is in fact easy to 
answer if ip is a Gaussian function because then its Wigner transform is pro- 
portional to a Gaussian e~~i S Sz ' z where S is a symplectic matrix uniquely 
determined by tp. It follows that there is a one-to-one correspondence be- 
tween Gaussians and the sets S T Sz ■ z < H. We have called these sets 

"quantum blobs" in [13 H3 [H E2 E3 El HB] ; the interest of these Quo- 
tum blobs comes from the fact that they represent minimum uncertainty 
sets in phase space. 

The Gaussian function 

tt„(s) = e~ x2 / 2h ; (1) 

is the (unnormalized) ground state of the one-dimensional harmonic oscil- 
lator with mass and frequency equal to one: HUfo = Eq^q where Eq = hh 
and 

s = \(-*&>^) < 2 > 

This operator is the quantization of the classical oscillator Hamiltonian 

H(x,p) = ^(p 2 + x 2 ) (3) 

The set Qq defined by the inequality H < Eq is the interior of the energy 
hypersurface H < Eq; it is the disk p 2 + x 2 < h with radius Rq = \f%. 
Let us now consider the iV-th excited state of the operator H; it is the 
(unnormalized) Hermite function 

q, N ( x ) = e - x2 l 2h H N {x/Vh) (4) 

where 

H N {x) = {-\) n e*£ r e-* (5) 

is the iV-th Hermite polynomial. It is a solution of H^n = (N + ^) H^n 
and the set £1n defined by the inequality H < En = (2N + 1)H is again a 

disk, but this time with radius = <J (N + |) h. 

In this paper we introduce a non-trivial extension of the notion of "quan- 
tum blob" we defined and studied in previous work. Quantum blobs are 
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deformations of the phase space ball \x\ 2 + \p\ 2 < h by translations and lin- 
ear canonical transformations. Their interest come from the fact that they 
provide us with a coarse-graining of phase space different from the usual 
coarse graining by cubes with volume ~ h n commonly used in statistical 
mechanics. They appear as space units of minimum uncertainty in one- 
to-one correspondence with the generalized coherent states familiar from 
quantum optics, and have allowed us to recover the exact ground states 
of generalized harmonic oscillators, as well as the semiclassical energy lev- 
els of quantum systems with completely integrable Hamiltonian function, 
and to explain them in terms of the topological notion of symplectic capac- 
ity [241 130] originating in Gromov's [23] non-squeezing theorem (alias "the 
principle of the symplectic camel"). Quantum blobs, do not, however, allow 
a characterization of excited states; for instance there is no obvious rela- 
tion between them and the Hermite functions. Why this does not work is 
easy to understand: quantum blobs correspond to the states saturating the 
S chr odinger-Rob ertson inequalities 



as is well-known [21] the quantum states for which all these inequalities 
become equalities are Gaussians, in this case precisely those who are them- 
selves the ground states of generalized harmonic oscillators. As soon as one 
consider the excited states the corresponding eigenfunctions are Hermite 
functions and for these the inequalities © are strict. The way out of this 
difficulty is to define new phase space objects, the "Fermi blobs" of the title 
of this paper. Such an approach should certainly be welcome in times where 
phase space is beginning to be taken seriously (see the recent review paper 



1.2 How We Will Do It 

We will show that a complete geometric picture of excited states can be 
given using an idea of the physicist Enrico Fermi in a largely forgotten pa- 
per [H] from 1930. Fermi associates to every quantum state \£ a certain 
hypersurface gp(x,p) = in phase space. The underlying idea is actually 
surprisingly simple. It consists in observing that any complex twice contin- 
uously differentiable function ^f(x) = R(x)e z ®( x " fi (R(x) > and <&(x) real) 
defined on M. n satisfies the partial differential equation 



(AXjfiAPj) 2 > A(Xj,Pj) 2 + \h 2 , 1 < j < n; 



(6) 
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where V 2 , is the Laplace operator in the variables xx,...,x n (it is assumed 
that R(x) ^ for x in some subset of M n ). Performing the gauge transfor- 
mation —iK\J x — > —iKV x — Vi;$, this equation is in fact equivalent to the 
trivial equation 

-H 2 V 2 x + h 2 ^-^R = 0. (8) 
The operator 

g$=(-ihV x -V x $) 2 + h 2 ^ (9) 

appearing in the left-hand side of Eqn. (J7]) is the quantisation (in every 
reasonable physical quantisation scheme) of the real observable 

9F(x,p) = (p-V x ^) 2 + h 2 ^- (10) 

and the equation g-p{x,p) = in general determines a hypersurface Tip in 
phase space M. 2 ™ p which Fermi ultimately identifies with the state ^ itself. 
The remarkable thing with this construction is that it shows that to an 
arbitrary function it associates a Hamiltonian function of the classical 
type 

H = (p - V^) 2 + V (11) 

even if ^ is the solution of another partial (or pseudo-differential) equation. 
We notice that when 1 F is an eigenstate of the operator = then 
gp = H — E and Tip is just the energy hypersurface H(x,p) = E. 

Of course, Fermi's analysis was very heuristic and its mathematical 
rigour borders the unacceptable (at least by modern standards). Fermi's 
paper has recently been rediscovered by Benenti [2] and Benenti and Strini 
[3], who study its relationship with the level sets of the Wigner transform 
of I'. 

Notation 1 The points in configuration and momentum space are written 
x = (xi, ■~,x n ) and p = (pi,--,p n ) respectively; in formulas x an p are 
viewed as column vectors. We will also use the collective notation z = (x,p) 

for the phase space variable. The matrix J = (^j (® an< ^ I the n x n 

zero and identity matrices ) defines the standard symplectic form on the phase 
space M. 2n via the formula o~(z,z') = Jz • z' = p ■ x' — p' ■ x. We write 
h = h/2n, h being Planck's constant. The symplectic group is denoted by 
Sp(2n,M): it is the multiplicative group of all real 2n x 2n matrices S such 
that a(Sz, Sz') = a(z, z') for all z, z' . 
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2 Symplectic Capacities and Quantum Blobs 



To generalize the discussion above to the multi-dimensional case we have to 
introduce some concepts from symplectic topology. For a review of these 
notions see de Gosson and Luef [22] , 

2.1 Symplectic Capacities 
Intrinsic symplectic capacities 

An intrinsic symplectic capacity assigns a non- negative number (or +oo) 
c(Q) to every subset Q, of phase space M 2n ; this assignment is subjected to 
the following properties: 

• Monotonicity: If ft C ft' then c(Q) < c(fi'); 

• Symplectic invariance: If / is a canonical transformation (linear, 
or not) then c(/(fi)) = c(£l); 

• Conformality: If A is a real number then c(XQ) = \ 2 c(Q); here XQ 
is the set of all points \z when z G f2; 

• Normalization: We have 



here B 2n (R) is the phase-space ball \x\ 2 + \p\ 2 < R 2 and Z 2n (R) the 



Let c be a symplectic capacity on the phase plane M. 2 . We have c(f2) = 
Area(f2) when Q is a connected and simply connected surface. In the general 
case there exist infinitely many intrinsic symplectic capacities, but they all 
agree on phase space ellipsoids as we will see below. The smallest symplectic 
capacity is denoted by c m \ n ("Gromov width"): by definition c m i n (0) is the 
supremum of all numbers ttR 2 such that there exists a canonical transfor- 
mation such that f(B 2n (R)) C O. The fact that c m ; n really is a symplec- 
tic capacity follows from a deep and difficult topological result, Gromov's 
|23j symplectic non-squeezing theorem, alias the principle of the symplectic 
camel. (For a discussion of Gromov's theorem from the point of view of 
Physics see de Gosson |17j . de Gosson and Luef [22].) Another useful exam- 
ple is provided by the Hofer-Zehnder |24j capacity c HZ . It has the property 



c(B 2n (R)) = irR 2 = c(Z 2n (R)); 



(12) 



phase-space cylinder x 2 + p 
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that it is given by the integral of the action form pdx = p\dx\ + • • • + p n dx n 
along a certain curve: 

c HZ (ft) = I pdx (13) 

■* Tmin 

when Q, is a compact convex set in phase space; here 7 m i n is the shortest 
(positively oriented) Hamiltonian periodic orbit carried by the boundary d£l 
of Q. This formula agrees with the usual notion of area in the case n = 1. 

It turns out that all intrinsic symplectic capacities agree on phase space 
ellipsoids, and are calculated as follows (see e.g. [16] [22], [23] ) . Let M be a 
2n x 2n positive-definite matrix M and consider the ellipsoid: 

n MjZ0 : M(z - z ) 2 < 1. (14) 

Then, for every intrinsic symplectic capacity c we have 

c(n M , Z0 ) = vr/A max (15) 

where A^ ax = is the largest symplectic eigenvalue of M. The symplectic 
eigenvalues of a positive definite matrix are defined as follows: the matrix 
JM (J the standard symplectic matrix) is equivalent to the antisymmetric 
matrix M 1 ' 2 JM 1 ' 2 hence its 2n eigenvalues are of the type izAf,.., ±iA^ 
where AJ > 0. The positive numbers AJ,.., A^ are called the symplectic 
eigenvalues of the matrix M. 

In particular, if X and Y are real symmetric n x n matrices, then the 
symplectic capacity of the ellipsoid 

n (AB) : Xx 2 +Y P 2 <1 (16) 

is given by 

C (^(A,B)) = 7r/ V Amax (17) 

where A max is the largest eigenvalue of AB. 
Extrinsic symplectic capacities 

The definition of an extrinsic symplectic capacity is similar to that of an 
intrinsic capacity, but one weakens the normalization condition f)12|) by only 
requiring that: 

• Nontriviality: c(B 2n (R)) < +00 and c(Z 2n (R)) < +00. 
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In [6] Ekeland and Hofer denned a sequence cf , cf H 



..,c 



EH 



trinsic symplectic capacities satisfying the nontriviality properties 



cf{B^{R)) 



k + n — 1 



ttR 



<f H ( 



Z 2n (i?)) = yfcvri? 2 



...of ex- 



(18) 



Of course cf H is an intrinsic capacity; in fact it coincides with the Hofer- 
Zehnder capacity on bounded convex sets f2. We have 



EH 



EH, 



(19) 



The Ekeland-Hofer capacities have the property that for each k there exists 
an integer N > and a closed characteristic 7 of d£l such that 



cf H (fl) = N 



pdx 



(20) 



(in other words, cF (fl) is a value of the action spectrum [5] of the bound- 
ary <9f2 of Q); this formula shows that c]? H (r2) is solely determined by <9f2; 
therefore the notation c^(dQ) is often used in the literature. The Ekeland- 
Hofer capacities c)? H allow us to classify phase-space ellipsoids. In fact, the 
non-decreasing sequence of numbers c^(0,m) is determined as follows for an 
ellipsoid £1 : Mz-z < 1 (M symmetric and positive-definite): let (Af, A^) 
be the symplectic eigenvalues of M; then 



: k = 1, 2, ...} = {iVvrA^ : j = 1, n; iV = 0, 1, 2, ...}. 



(21) 



Equivalently, the increasing sequence c^ H (f2) < C2 H (0) < ■ • • is obtained by 
writing the numbers in increasing order with repetitions if a number 

occurs more than once. 



2.2 Quantum Blobs 

By definition a quantum blob QB 2n (zo, S) is the image of the phase space 
ball B 2n (S~ 1 ZQ, y/h) : \z — S~ 1 zq\ < yH by a linear canonical transformation 
(identified with a symplectic matrix 5). A quantum blob is thus a phase 
space ellipsoid with symplectic capacity ttH = ^h, but it is not true that, 
conversely, an arbitrary phase space ellipsoid with symplectic capacity is 
a quantum blob. One can however show (de Gosson [HI [13 [16], de Gosson 
and Luef [22 j) that such an ellipsoid contains a unique quantum blob. One 
proves (ibid.) that a quantum blob QB 2n (zo, S) is characterized by the two 
following equivalent properties: 



7 



The intersection of the ellipsoid QB 2n (zo, S) with a plane passing 
through zq and parallel to any of the plane of canonically conjugate 
coordinates Xj,pj in M. 2n is an ellipse with area ^h; 



• The supremum of the set of all numbers irR 2 such that the ball B 2n (\/ r R) 
\z\ < R can be embedded into QB 2n (zo, S) using canonical transfor- 
mations (linear, or not) is \h. Hence no phase space ball with radius 
R > Vh can be "squeezed" inside QB 2n (zo, S) using only canonical 
transformations. 

It turns out (de Gosson |16j) that in the first of these conditions one 
can replace the plane of conjugate coordinates with any symplectic plane (a 
symplectic plane is a two-dimensional subspace of M^™ 011 which the restric- 
tion of the symplectic form a is again a symplectic form) . There is a natural 
action 

Sp(2n,R) x QB(2n, R) — ► QB(2n,R) 
of the symplectic group on quantum blobs. 

3 Generalized Coherent States 

3.1 The Fermi Function of a Gaussian 

We next consider arbitrary (normalized) generalized coherent states 



1 \ n/4 

^x,y(x)=[—\ (detX)V*exp 



--^(X + iY)x ■ x 



(22) 



where X and Y are real symmetric nxn matrices, and X is positive definite. 
Setting $(x) = —\Yx ■ x and R(x) = exp (— ^Xx ■ x) we have 

= - y * • ^ = >* + ?A, (23) 
hence the Fermi function of ^x,Y is the quadratic form 

g F (x,p) = (p + Yx) 2 + X 2 x ■ x - HTvX. (24) 
We can rewrite this formula as 

g F (x,p) = M F z ■ z - HTvX (25) 



S 



(z = (x,p)) where Mp is the symmetric matrix 



Mf=[ x 2 + y* y\ m 



Y I 

A straightforward calculation shows that we have the factorization 

M F = S T ( X Q ° x "jS (27) 
where S is the symplectic matrix 

S =(x-i/2 Y X V)- (28) 

It turns out -and this is really a striking fact!- that Mp is closely related 
to the Wigner transform 

W* x , Y (z) = e-*™*x,Y{x + \yWxA x ~ b) d y ( 29 ) 

of the state ^x,Y because we have 

W* x ,y(z) = (J^j exp (~Gz ■ z^j (30) 
where G is the symplectic matrix 

G = S T S=( X+ /_f Y lY ^) (31) 

(see e.g. |16[ I27|). When n = 1 and ^x,y( x ) = ^o( x ) the fiducial coherent 
state (HJ) we have S~ 1 D~ 1 / 2 S = I and Tr X = 1 hence the formula 

/ 1 \ 1/4 1 r i 

W^ (z) = — - exp ——Mpz ■ z 
\7rnJ e [ n 

already observed by Benenti and Strini in [3]. 
3.2 Geometric Interpretation 

Recall (formula (fT5|) ) that the symplectic capacity c(f2) of an ellipsoid Mz ■ 
z < 1 (M a symmetric positive definite 2n x 2n matrix) is given by 

c(fi) = ir/\^ x (32) 
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where A max = max{Af,.., A^}, the AJ being the symplectic eigenvalues of 
M. We denote by fl F the phase space ellipsoid defined by g-p(x,p) < 0, that 
is: 

n F : M F z • z < hTrX; 

it is the ellipsoid bounded by the Fermi hypersurface T~L F corresponding to 
the generalized coherent state fjy, Let us perform the symplectic change 
of variables z' = Sz; in the new coordinates the ellipsoid tt F is represented 
by the inequality 

Xx' ■ x' + Xp' ■ p' < hTr X (33) 

hence c(Q F ) equals the symplectic capacity of the ellipsoid (J33j) . Applying 
the rule above we thus have to find the symplectic eigenvalues of the block- 

'X 0' 



diagonal matrix I ^ ; a straightforward calculation shows that these 

are just the eigenvalues wi, ui n of X and hence 

c(n F ) =vr?iTrX/a; max (34) 
where w max = maxjwi, ...,a; n }. In view of the trivial inequality 

n 

Umax < TrX = y~]u)j < nAw max (35) 
i=i 

we have 

\h < c(Q F ) < ^. (36) 

An immediate consequence of the inequality \h < c(J2f) is that the Fermi 
ellipsoid Op of a generalized coherent state always contains a quantum blob; 
this is of course consistent with the uncertainty principle. 

Notice that when all the eigenvalues ujj are equal to a number to then 
e(f2p) = nh/2; in particular when n = 1 we have c(f2p) = h/2 which 
is exactly the action calculated along the trajectory corresponding to the 
ground state. This observation leads us to the following question: what is 
the precise geometric meaning of formula (|34p ? Let us come back to the 
interpretation of the ellipsoid defined by the inequality (f33|h We have seen 



that the symplectic eigenvalues of the matrix I ^ ^.1 are precisely the 



X y 

eigenvalues Uj, 1 < j < n, of the positive-definite matrix X. It follows that 
there exist linear symplectic coordinates (x",p") in which the equation of 
the ellipsoid fip takes the normal form 

n n 
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whose quantum-mechanical interpretation is clear: dividing both sides by 
two we get the energy shell of the anisotropic harmonic oscillator in its 
ground state. Consider now the planes V\, V 2 , ■■,V n of conjugate coordinates 
(x\,pi), (x 2 ,p 2 ) v .., ( x n,Pn)- The intersection of the ellipsoid ftp with these 
planes are the circles 



Ci luxixf +p'{ 2 ) < ^hujj 



3=1 
n 



C 2 : u 2 {x' 2 a + p'f) < ^huj 
3=1 



C n : U n (x'n +Pn ) < ^2 t™3' 
3=1 

Formula (|34p says that c(J1f) is the area of the circle Cj with smallest radius, 
and this corresponds to the index j such that Uj = o; max . This is of course 
perfectly in accordance with the definition of the Hofer-Zehnder capacity 
c HZ (f2p) since all symplectic capacities agree on ellipsoids. We are now 
led to another question: is there any way to describe topologically Fermi's 
ellipsoid in such a way that the areas of every circle Cj becomes apparent? 
The problem with the standard capacity of an ellipsoid is that it only "sees" 
the smallest cut of that ellipsoid by a plane of conjugate coordinate. The 
way out of this difficult lies in the use of the Ekeland-Hofer capacities c™ 
discussed above. To illustrate the idea, let us first consider the case n = 2; 
it is no restriction to assume cj% < uj 2 - If oj\ = oj 2 then the ellipsoid 

wi Of + p'( 2 ) + u 2 {xf + pf )<hw 1 + nw 2 (38) 



is just the ball B 2 (y2h) whose symplectic capacity is 2irh = h. Suppose 
now that oj\ < uj 2 . Then the Ekeland-Hofer capacities are the numbers 

ttK, , ffli, > 2irh 2ttH 

— (cji + w 2 ), — (wi + u 2 ), (wi + w 2 ), (wi + u 2 ), .... (39) 

U3 2 UJ\ U3 2 D\ 



and hence 



c f H (fi F ) = c(n F ) = — ( Ul + W 2 ). 



0J 2 

What about cf H (f2p)? A first glance at the sequence ([39]) suggests that we 
have 

c f H (r2 F ) = — ( W i +W2 ) 
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but this is only true if uj\ < UJ2 < 2cji because if 2cji < ^2 then (cji + 
0O2)/(jO2 < (oji +uj2)/oji so that in this case 

The Ekeland-Hofer capacities thus allow a geometrical classification of the 
eigenstates. 

4 Fermi Function and Excited States 

The generalized coherent states can be viewed as the ground states of a 
generalized harmonic oscillator, with Hamiltonian function a homogeneous 
quadratic polynomial in the position and momentum coordinates: 

H(x,p) = ^aijPiPj + bijPiXj + CijXiXj. 

Such a function can always be put in the form 

H(z) = ^Mz ■ z (40) 

where M is a symmetric matrix (the Hessian matrix, i.e. the matrix of 
second derivatives, of H). We will assume for simplicity that M is positive- 
definite; we can then always bring it into the normal form 

n 

*(*) = Ef(*?+*S) 

3=1 

using a linear symplectic transformation of the coordinates (symplectic di- 
agonalization) : there exists a symplectic matrix S (depending on M) such 
that 

S T MS = D=(^ °) (41) 

where SI is a diagonal matrix whose diagonal entries consist of the symplectic 
spectrum u>i, ...,u n of M. Thus, we have K{z) = H(Sz), or, equivalently, 

H(z) = KiS^z) (42) 

The ground state of each one-dimensional quantum oscillator 
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is the solution of Kj*$> = ^hioj^, it is thus the one-dimensional fiducial coher- 
ent state (7rfr) _1 / 4 e -x ' 2//2?l . It follows that the ground state of K = ^ . Kj 

is the tensor product of n such states, that is ^!$(x) = (7r/i)~ n//4 e _ ' z ' 2//2fi , the 
fiducial coherent state (H|). Returning to the initial Hamiltonian H we note 
that the corresponding Weyl quantisation H satisfies, in view of Eqn. (142ft 
the symplectic covariance formula H = SKS^wheie S is any of the two 
metaplectic operators corresponding to the symplectic matrix S (see the 
Appendix). It follows that the ground state of H is given by the formula 
* = 5*o- 

The case of the excited states is treated similarly. The solutions of the 
one-dimensional eigenfunction problem Kj ^ = are given by the Hermite 
functions 

m N {x) = e - x2 / 2h H N (x/Vh) (43) 

with corresponding eigenvalues En = (N + ^)hujj. It follows that the solu- 
tions of the n-dimensional problem = E^f are the tensor products 

tf (JV) = ®^N 2 <8> • • -® #iV n (44) 

where (N) = (N\,N2, ...,N n ) is a sequence of non- negative integers, and the 
corresponding energy level is 

n 

E (N) = + h)^r (45) 

i=i 

This allows us to give a geometric description of all eigenfunctions of the 
generalized harmonic oscillator, corresponding to a quadratic Hamiltonian 
(jlOj) . We claim that: 

Let be an eigenfunction of the operator 

H = (x, -ihV x )M(x, -ihV x ) T . (46) 

The symplectic capacity of the corresponding Fermi blob Qp is 

n 

c(n F ) = J2(N 3 + h)h (47) 
i=i 

where the numbers N\, N2, N n are the non-negative integers 
corresponding to the state (1441) of the diagonalized operator K = 
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APPENDIX: The Metaplectic Group 



The symplectic group Sp(2ra, R) has a covering group of order two, the meta- 
plectic group Mp(2n,R). That group consists of unitary operators (the 
metaplectic operators) acting on L 2 (IR n ). There are several equivalent ways 
to describe the metaplectic operators. For our purposes the most tractable 
is the following: assume that S £ Sp(2n,R) has the block-matrix form 

S = ^\ with det B^O. (Al) 

The condition det B ^ is not very restrictive, because one shows (de Gos- 
son JD1 ES ES], Littlejohn [27]) that every S G Sp(2ra,R) can be written 
(non uniquely) as the product of two symplectic matrices of the type above; 
moreover the symplectic matrices arising as Jacobian matrices of Hamilto- 
nian flows determined by physical Hamiltonians of the type "kinetic energy 
plus potential" are of this type for almost every time t. To the matrix (|A1|) 
we associate the following quantities (de Gosson fL0\ [16]): 

• A quadratic form 

W(x, x') = -DB^x ■ x - B~ l x ■ x 1 + -B~ x Ax' ■ x'; (A2) 

the matrices DB^ 1 and B~ 1 A are symmetric because S is symplectic; 

• The complex number A(PI^) = i m \J\ det B~ 1 \ where m ("Maslov in- 
dex") is chosen in the following way: m = or 2 if det B' 1 > and 
m = 1 or 3 if det B' 1 < 0. 

The two metaplectic operators associated to S are then given by 

= (2^l) n/2 A(W0 J e* w( - x > x '^{x')d n x'. (A3) 

The fact that we have two possible choices for the Maslov index is directly 
related the fact that Mp(2n, R) is a two-fold covering group of the symplectic 
group Sp(2n,R) pH EDI EH i • 

The main interest of the metaplectic group in quantization questions 
comes from the two following (related) "symplectic covariance" properties: 

• Let f be a square integrable function (or, more generally, a tempered 
distribution), and S a symplectic matrix. We have 

W^(S- l z) = W(SV)(z) (A4) 

where S is any of the two metaplectic operators corresponding to S; 
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• Let H be the Weyl quantisation of the symbol (= observable) H. Let 
S be a symplectic matrix Then the quantisation of K{z) = H(Sz) is 
K = S~ l HS where S is again defined as above. 
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